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Abstract 

We study the critical effect of an intermittent social distancing strategy on the propagation 
of epidemics in adaptive complex networks. We characterize the effect of our strategy in the 
framework of the susceptible-infected-recovered model. In our model, based on local information, 
a susceptible individual interrupts the contact with an infected individual with a probability a 
^ ■ and restores it after a fixed time ti,. We find that, depending on the network topology, in our 

i-^ . social distancing strategy there exists a cutoff threshold cTc beyond which the epidemic phase 

disappears. Our results are supported by a theoretical framework and extensive simulations 
of the model. Furthermore we show that this strategy is very efficient because it leads to a 
"susceptible herd behavior" that protects a large fraction of susceptibles individuals. We explain 
our results using percolation arguments. 

PACS numbers: 89.75.-k, 64.60.aq, 64.60.ah 
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I. INTRODUCTION 



The study of the topology of complex networks, and the dynamical processes that 
use these networks as substrate to spread, has recently generated great interest in the 
scientific community jl-^. In the past, studies of dynamic processes such as the spread- 
ing of rumors, opinions, and diseases on static networks, were concentrated on how their 
topology affects these processes jl-3|. However, it is known that these processes evolve 
on top of networks where the topology changes with time j4|. As a consequence, recently 
many researchers began to study dynamic networks, and the interplay between the dy- 
namic process and the network dynamics. Those networks in which the topology changes 
regardless of the process taking place on top of them are called evolutive networks, while 
those networks that change their topology to mitigate or promote these processes are 
called adaptive networks 0, Adaptive networks have been investigated in many dis- 
ciplines such as social sciences, epidemiology, biology, etc. ji, ^ 7|. In these networks 
there is a coevolution between the link dynamics and the state of the nodes which leads 
to a collective phenomenon on adapting networks. As an example, in a network of routes, 
where the nodes are cities and the links are the routes connecting them, some overloaded 
paths become dysfunctional and new paths between cities are built to avoid the conges- 
tion. In the analysis of opinion formation on social networks, nodes usually tend to rewire 
or break their links with individuals with different opinions, leading in some models to a 
network fragmentation into components or clusters in which all members have the same 
opinion j7-9|. In biological networks, such as the vascular system, after an arterial occlu- 
sion collateral vessels grows in order to increase the blood flow to neighboring tissues jlO|. 
Similarly, in the widely studied epidemic models on static networks, adaptive pro cesses 
are used to model strategies that reduce the impact of the disease spreading jsl. Illi 12 1. 

One of the most popular models in epidemiology that reproduces seasonal diseases is 
the susceptible-infected-recovered model (SIR) jll.1l3|. where individuals can be in one of 
three states, susceptible (S), infected (J) or recovered (i?). In its classical formulation, 
an infected node infects a susceptible neighbor with probability /3 and recovers with a 
certain fixed probability, which implies an exponential distribution of times for which 
individuals remain infected. However, this distribution is rarel y r ealistic and for most 
seasonal diseases it has a sharp peak around an average value As a consequence. 
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some studies 
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17| have used a different version of the SIR model in which an infected 



individual recovers after a fixed time t^, called the recovery time. 

It is well known that in the SIR model on static networks, the size of the infection 
is governed by the effective probability of infection or transmissibility T of the disease, 
where T = 1 — (1 — /?)*''. In turn, it was shown that this model can be mapped into 



a link percolation process 
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18|, where T plays the role of the link occupancy proba- 
bility p in percolation. In a percolation process, there is a critical probability where 
the finite cluster size distribution behaves as ~ in the thermodynamic limit. 
Above this threshold a "giant component" appears. As a consequence of the mapping 
between percolation and the SIR model; in the latter there is an epidemic threshold at 
Tc = Pc below which the disease is an outbreak were the infection reaches a small frac- 
tion of the population, which is equivalent to having finite clusters in percolation, while 
above Tc an epidemic develops corresponding to the emergence of a percolating giant 



component 



19 



20| . This threshold, in uncorrelated static networks, depends only on 



the degree distribution P{k), where k is the degree or the number of links that a node 
can have. In particular, for Erdos-Renyi (ER) networks, P{k) = e"'^^'' {k)^ /k\ where {k) 
is the mean connectivity, the threshold is Tc = l/{k). However, in pure scale- free (SF) 
networks P{k) ~ k^^, where A is the broadness of the distribution, in the thermodynamic 
limit Tc — )■ for A < 3, which means that the epidemic spreads for any value of T. This 
last result indicates that highly heterogeneous networks, such as many theoretical social 
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networks, are very likely to develop an epidemic 

While it is well known how the topology affects the SIR process in static networks, 
there is very little literature about this model on adaptive networks. Recently, Lagorio 
et al. Q studied two different strategies to mitigate the spread of a disease in an adaptive 
SIR model. In that model a susceptible node disconnects the link with an infected individ- 
ual with a rewiring probability w, and creates a link with another susceptible node. The 
authors found that there is a phase transition at a critical rewiring threshold Wc separating 
an epidemic from a non-epidemic phase, which can be related to static link percolation. 
A feature of this rewiring process is that links between susceptible individuals can be 
established independently of their previous relationship. In that strategy, the nodes have 
no memory because two individuals can be connected independently of their past. Even 
though this adaptive process could be representative of casual contacts between individu- 
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als such as those generated in pubhc buildings hke shopping, theaters, etc., that strategy 
will not work for other types of interaction such as friendship and working partners, where 
individuals preserve their closer contacts. Therefore, if an individual is separated from its 
closer neighbors, it will tend to reconnect with them, at some time, more often than with 



an unknown individual. On the other hand, recently Wang et al. [22| and Van Segbroeck 
et al. {231 proposed strategies to stop the spread in the SIR model where the susceptible 
individuals driven by fear disconnect their links with their neighbors, infected or not, 
without creating a new link. This strategy reduces the number of contacts permanently, 
which is efficient but also very inconvenient from an economical point of view. 

In this paper we propose a strategy based on "intermittent social distancing" in adaptive 
networks and study its efficiency in stopping the spread of diseases in theoretical and 
real networks. We found theoretically that in our model there exists a cutoff threshold 
that prevents an epidemic phase. Our results are supported by extensive simulations. 
Moreover, we found that the intermittent social distancing strategy is efficient to protect 
a large susceptible cluster. The paper is organized as following: in Sec. [TTl we derive 
the theoretical transmissibility for our model and show how our strategy diminishes the 
epidemic phase. In Sec. IIIH we show how the epidemic size is reduced with our strategy 
and the agreement between our theoretical approach and the simulations. In Sec. IIVI we 
present a study of the "susceptible herd behavior" that we use as a criterion to evaluate 
the effectiveness of our strategy. In Sec. El we present our conclusions. 



II. ANALYTICAL APPROACH 



We propose a SIR model in an adaptive network where an infected individual transmits 
the disease to a susceptible neighbor with probability (3 and, if he fails, with probability 
cr the susceptible individual breaks the link with the infected one for a period tb. Thus, 
the effective probability of breaking a link is (1 — /3)cr. After a time th both nodes are 
reconnected and the process is repeated until the infected node recovers at a fixed time 
tj. > tfe, i.e., there is an intermittent connection between a susceptible node and its 
infected neighbor. This mimics a behavioral adaptation of the society to avoid contacts 
with infected individuals by imposing a social distancing during one or more periods 
of duration ti,. Notice that in our model a susceptible node breaks its links using only 
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local information and not global knowledge as provided by communication media. In our 
model, the time is increased by 1 after every infected node tries to infect its neighbors and 
the updates are done after each time step. In this process the dynamic transmissibility 
T(/3, 0", tr, tb) = Tcr can be written as 

r, = ^/3(l-/3)""i(l-ar"i + /3 J2 ^in,h,a,(3y, (1) 

n=l n=ti,+2 

where 0(n, U, a, (3) is given by 

0(n,t,,a,/3)= J2 (''"''^'"^)^"(l-^r^~"^*^^^Hl-/3r'""*^ (2) 

n=l V ^ / 

and [■ ■ ■ ] denotes the integer part function. The first term of Eq. ([1]) is the probability that 
a node in the I state transmits the disease to any neighbor node in the S state (before it 
recovers), considering that the link S-I has never been broken. The second term represents 
the probability for a node in the / state to transmit the disease to an S node after that pair 
has been disconnected u times for a period t^. The binomial coefficient of the second term 
takes into account the number of ways to arrange u intermittent disconnected periods 
before the susceptible becomes infected individual at time n. Notice that uth is the total 
times that the pair S-I is broken and represent a temporal social distancing. In Table [U 
we illustrate the element n = 8 of the second term of Eq. ([1]) with tr = 10 and th = 2. 



5 



TABLE I: Disconnected periods for a pair S-I with tr = 10 (recovery time), th = 2 (dis- 
connection period), and n = 8 (time of infection). The first column represents the number 
of disconnected periods u before n = 8, the second column is a typical configuration, the 
third column is the probability of that configuration, and the fourth column is the number 
of ways to arrange u disconnected periods. In the second column, each cell corresponds to 
a unit time. The white cells represent the time unit where a link between the S and the / 
node exists, the gray ones correspond to the disconnection period, and in the black cells 
there is no dynamic for the pair S-I because the S has been infected and now the pair 
becomes I-I. Notice that initially the link cannot be broken because this disconnection 
happens only after the I individual fails to infect the susceptible one, with probability 
(1 — /3). Similarly, two disconnection periods must be separated by at least a white cell. 
When the infection occurs at time n, the maximum number of disconnected periods is 
u = [{n-l)/{h+l)]. 

u Example Probability Binomial Coefficient 

t, 

u = l (3 a{l-af{l- I3f (8-2-I) = 5 

^= 1 



It is known that the disease becomes an epidemic if the basic reproductive number 
-Ro > 1; where Rq is the number of secondary infections. For uncorrelated networks, Rq 
is related to the connectivity distribution P{k) through the branching factor k, 

R^ = {k- 1)T,, (3) 

where k = {k'^)/{k), and (k) and {k"^) are the first and the second moments of P{k) 
respectively. The branching factor is a measure of the heterogeneity of the network that 



diverges for SF networks with A < 3 in the thermodynamic limit because {k"^) — 00 2l|. 
Then the critical cutoff (Xc, above which the disease dies, as a function of /3, tr, and th, 
can be found through the condition Rq = 1, which yields 

Ta. = = Tc, (4) 
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2ll . |24| where is the critical transmissibihty for the SIR model in static networks, and 
consequently our dynamic process in the steady state is related to a static topological 
property of the network. This is expected due to the fact that in our model the network 
topology does not change globally in the characteristic time scale of the disease spreading, 
and hence our strategy can be understood as an SIR model on a static network but with 
a transmissibihty Tc = T^^. 

From Eq. ([1]) it is straightforward that in the limit a— j-O, To- = T= 1 — (1 — /?)*''. On 
the other hand, when a — t- 1, as the only terms that survive in Eq. are those which 
fulfills the condition n — 1 — u{ti, + 1) = 0, we obtain 



T 

■J- rr 



(5) 

In Fig.lll we plot the plane a—T [T = T{(j = 0)] in order to show that with our strategy 
the epidemic phase is reduced compared to the static case. Notice that for tb = tr/2 the 
epidemic phase shrinks substantially compared to the case = 1. 
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FIG. 1: Plot of the epidemic phases in the plane a — T for = 20 and static = 0.25, 
where T corresponds to the transmissibility in a non adaptive network . The dashed 
lines correspond to the critical dynamic transmissibility To-^ for (from left to right) t;, = 1, 
tb = tr/2 and ti, = tr — l. For tfe = 1 and o" = 1, the maximum transmissibility for the static 
SIR model for which the epidemic phase is T = 1 - (1 - T^yr/[{tr+i)/2] or T ^ 1 - {l-Tcf, 
and the tr dependence disappears. 

In the figure, the light-gray area, delimited between the curves corresponding to the 
blocking periods tt = 1 and tb = tr — 1, displays the region of parameters controlled by 
the intervention strategy. We can see that even the milder intervention = 1 expands 
the epidemic-free area compared to the static case. On the other hand, the stronger 
intervention t/, = tr — 1 drastically shrinks the epidemic phase as a increases. From the 
figure, we can also see that for Tc ^ 0, as the social distancing increases the epidemic-free 
phase increases; meanwhile, for a theoretical SF network with A < 3, as Tc = 0, when 

— 7- oo there is no epidemic-free phase for any values of a and tb- However, real networks 
are finite and they are not generally pure uncorrelated SF networks which implies that 
our intermittent social distancing strategy could be applied in these networks. 

In order to determine how the heterogeneity of the network affects the strategy per- 
formance, in Fig. [2] we plot a phase diagram in the plane tb — a for fixed tr and different 
values of k. 
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FIG. 2: Phase diagram for a and tb for /3 = 0.05 and = 20 for different heterogeneities 
K. The dashed hnes with circles represent the interface between the non-epidemic (right) 
and epidemic (left) phases for the different values of k. Notice that Oc = 0.33 and 
Oc = 0.86 correspond to the critical cutoff probabilities where in some region the interface 
is a vertical line for k = 5 and k = 10, respectively. 



We can see that as the heterogeneity k increases, the social distancing tt and cr have 
to increase in order to prevent the epidemic phase. Surprisingly, in high heterogeneous 
networks, we find that for th > tr/2, the critical cutoff probability o"c = (5c is almost 
constant (see the Appendix with 

(7) 

i.e., for very heterogeneous networks above ti, = tr/2, does not change the transmis- 
sibility; then the best and least expensive strategy is to reconnect at tb = tr/2. This 
means that if we know the average duration of a disease tr, individuals can return to 
their activities with low risk just after half of the characteristic time. From Eq. ([7j), it is 
straightforward that for Oc = 1 there exists an upper value of k, = k'™ that depends only 



P 

(1 - B) 



K-2 + ^/{k - 2){k- 1) 
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on f3, where is given by 



jin, 2 - (1 - By 



where this strategy can be apphed. From Eq. the hmit = 1 can also be expressed 
in terms of the minimum critical transmissibility, 

= 2/3 - p\ (9) 

As a consequence, for very low infection rates, our strategy predicts an epidemic-free phase 
even for highly heterogeneous networks with a finite epidemic threshold. Notice that since 
real networks have degree correlations and clustering, then the relation Tc = — 1) 
does not hold and the only magnitude that matters is T, which has to be measured by 



the peak of the second giant component 
3q. ([9]) instead of Eq. ([8]). We found that for the condensed matter coauthorship network 



25| . In those cases, for cTc = 1 we have to use 
[or the condensed matter coauthorship network 

n 

26] (Tc = 0.026) and mathematics coauthorship network [27| (T^ = 0.050) the epidemic 
spread can be stopped in diseases with /3 < 0.013 and /3 < 0.025, respectively. Similar 
values of (3 were used for real networks by Kitsak et al. 



III. NUMERICAL RESULTS 



In our model at the initial stage, all nodes are susceptible and we infect a randomly 
chosen node (patient zero) in the biggest component of our network. From the patient zero 
the disease spreads to its neighbors according to the rules of our model described above: 
an infected individual transmits the disease to a susceptible neighbor with probability 
/3 and, if it fails, with probability a the susceptible individual breaks the link with the 
infected one for a period After a time both nodes are reconnected and the process 
is repeated until the infected node recovers at a fixed time tr > tb- The time is increased 
by 1 after every infected node tries to infect its neighbors and the updates are done after 
each time step. 

All our results are presented for tr = 20 but qualitatively all the results are the same 
for tr > 1. 
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FIG. 3: Mj{a,tb)/Mi{a = 0) vs. a in an ER network with (A;) = 4, iV = 10^ U = 20, for 
different values of U: t;, = 3 (A), tf, = 5 (o ), 4 = 10 (O), U = 15 (□) and U = 19 (Q) 
for P = 0.05 with original transmissibility T = 0.64 (a) and /3 = 0.025 with T = 0.40 
(b); Mj{a,tf,)/Mj{a = 0) vs. a in a SF network with A = 3.5 and minimal connectivity 
kmin = 2 for tr = 20, /3 = 0.075 with original transmissibility T = 0.79(c) and /3 = 0.05 
with original transmissibility T = 0.64 (d). In the insets we show an enlargement of 
the main plot, rescaled in the abscissa by the factor cTc, obtained from Eq. ([1]). Our 
simulations were averaged over 10^ realizations. 

In Fig. [3l we plot the relative epidemic size M//M/(0) = Mj{a;tb) /Mj{a = 0) as a 
function of a for ER and SF networks for different values of th and = 20. From the plot 
we can see that M/(cr) decreases as a and tb increase compared to the static case M/(0). 
We can also see that a critical probability exists, which can be obtained theoretically 
from Eq. (jlj), above which the disease dies. Then, depending on how virulent is the disease 
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without intervention, we can control in order to stop the spread. In the insets of the 
figures, we collapse all the curves using the value of cTc obtained from Eq. (jl]). The collapse 
close to the critical value shows the excellent agreement between the theoretical value of 
(Tc and the simulations. 



IV. SUSCEPTIBLE HERD BEHAVIOR 



With our intermittent social distancing strategy, the susceptible nodes dynamically 
reduce their contact with the nodes in the infected cluster, mitigating the spread of the 
disease, i.e., our strategy produces a resistance to the disease which we call "susceptible 
herd behavior". As a result of the coevolutive process, at the end of the spreading there 
is only one cluster composed of recovered individuals and, depending on one or more 
susceptible clusters that give raise to a cluster size distribution of susceptible individuals or 



"voids" 29|. In our model, the cluster size distribution of voids or susceptible individuals 
is important since the formation of a susceptible herd behavior induced by the network 
dynamics also measures how effective our strategy is to preserve a whole part of the 
society safe from the disease. Next, we derive the value of the transmissibility for which 
a susceptible crowding develops, i.e., the value below which our strategy is efficient. 



Describing the growth of an epidemic cluster as a Leath process [30|, |31| for a value 
of the link occupancy probability p = T„ and denoting by the probability that a 

cluster reaches the nth generation following a link, then the probability /oo(p) that a link 
leads to a giant component when — oo is given by 

/-(p) = E^[i-p/-(i>)f-\ (10) 

k=i ^ ' 

where fooip) is the solution of 

fM = ^-Gi{i-v fUv)). (11) 

and Gi(x) = Er=i kP{k)/{k)x''-\ 

When the "epidemic" cluster grows, the size of the void clusters is reduced as in a node 
dilution process, since when a link is occupied a void cluster loses a node and all its edges. 
Then foo{p) is the probability that a void cluster loses a node. If we denote by 1 — p'' the 
fraction of void nodes removed the following relation holds 32|, 

l-p^ = /oc(p). (12) 
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For node percolation, it is known that [33[ P^{p") + n1 = p" , where P^{p^) is the 
fraction of nodes in the giant void component and the number of finite void clusters of 
size s. The fraction of remaining nodes, below which the giant void cluster is destroyed 
corresponds to the critical probability of node void percolation = Pc- At this value 
^ooiPc) — 0; then '^s = p"^. This means that at the void transition, only a fraction p"^ 
of the nodes belong to void clusters. As a consequence, the fraction of links p* needed to 
reach this point fulfills 

pI = i-Up*). (13) 

Therefore, from Eqs. (fTT]) and fll3p 



G,{l-p\l-pl)), 



(14) 



where pi = Gi ((G''^)~^(l)) and p* is the solution of Eq. ( I14p . Notice that at the void 
transition G'i(l — p* + p^p*) = 1. A similar result was also found by Newman for a process 
of spreading of two pathogens j34|, employing a different approach. Using the mapping 
between our adaptive SIR model and link percolation, the transmissibility T* needed to 
create a giant component of crowded susceptible individuals is T* = p*. Thus, for a 
disease spreading in an ER network, the dynamical critical transmissibility for the giant 
susceptible cluster is given by T* = —T„^ ln(To-^)/(l — T^J. 

In Fig. m we compare the cluster size distribution of susceptible individuals or voids 
nl for different values of the transmissibility where we distinguish outbreaks from epi- 



demics 
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16 



, |35| on an ER network with (k) = 4. 
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FIG. 4: Log-log plot of cluster size distribution of susceptible individuals for an ER 
network with (A;) = 4 {T* = 0.46) and N = 10^ at: = 0.64 (A), = 0.46 (Q) and 
To- = 0.40 (□). The dotted lines are a guide to the eye and the dashed line is the result 
of a power law fitting for T*. Notice that we add in the giant component. 



For To- > T* we obtain only small clusters of susceptible nodes that decay faster than 
exponentially. On the other hand, for T^- < T^ a giant susceptible cluster appears, meaning 
that our strategy is efficient producing large connected clusters of susceptible nodes that 
crowd in order to protect themselves. When with our strategy the transmissibility is 
reduced to T^- = T* the epidemic spreading slows down and the distribution decays as 
a power law with the same exponent r = 5/2 as in a void node percolation transition, 
in contrast with the distribution of infected sizes of outbreaks which at criticality goes 
as rig ~ s""^"*"^ as in a Leath process 30|. These results confirms the importance of the 
interplay between void nodes and link percolation in our model. . 



V. CONCLUSIONS 



In this paper we present a novel adaptive strategy based on intermittent social distanc- 
ing in an adaptive SIR model. In the intermittent social distancing strategy a susceptible 
individual breaks the link with the infected neighbor with a cutoff probability cr and then 
both individuals are reconnected after a time tb before the infected individual recovers at 
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tj.. Using the framework of percolation theory, we derive the dynamical transmissibility, 
and we find that there exists a critical cutoff Cc where the epidemic spread is stopped for 
non-highly heterogeneous networks. We show that in some real networks our intermittent 
social distance strategy could stop the epidemic spreading for not very virulent diseases. 
We find an excellent agreement between the theory and the simulations. For heteroge- 
neous networks, we find that a very high value of tb does not lead to a decrease of the 
epidemic, which implies that the less expensive strategy is to chose tb = tr/2. Finally 
we verify that our strategy reduces the transmissibility below a value where a susceptible 
crowding is produced. Any clever strategy used to reduce the disease from spreading 
should protect the population at the least economic cost. We believe that our strategy, 
which allows us to control the disconnection of periods through tb, is a very convenient 
strategy because it creates a susceptible herd cluster. Our present findings could be used 
as a support and reference guidance for the development of further strategies to stop 
diseases in real networks. 
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Appendix A 

For tb > [{tr — l)/2] ^ tr/2 and tb < tr — 1, in Eq. ([2]) the summation has only one term, 
because for tb > tr/2 we have only one period of disconnection, i.e., [{n — 1)/ {tb + 1)] = 1. 
Intuitively, if > tr/2 then in order to compute the transmissibility, we have to consider 
at most only one break period since, otherwise the transmission of the disease is not 
possible. Denoting by 5 = (1 — cr)(l — f3), then Eq. ([T]) reduces to 

= ^/3 5"-i + /3(l-/3)a J2 {n-h-m' 



rn— 2— tf. 



n=l n=ti,+2 

1-6 ^ (1 - a) (1 - 5)2 



[1 + 5^-4-1 ^S{tr -tb-l)+tb- tr)] . (Al) 
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Neglecting the higher powers of 6 (which hold for a — ?> 1, /3 — 1, or high values of tr —tf,), 
we obtain 

1-5+ (1-5)2 ■ ^ ^ 

Notice that T^- loses all the dependence on tr and as shown in Fig. |2j With this 
approximation and using the fact that = = 1/(1 — k), the critical cutoff ac = 5^c is 
given by 




^c = T^h^-l + A • (A3) 



If Eq. (HD holds, then 



/3 



1-/3 



K - 2 + v/(/t- l)(ft;-2) . (A4) 
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